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Abstract 

The Standard Model predicts a large width difference in the B s system. New 
physics can contribute significantly to the mass difference. If this contribution 
is CP violating, this will always result in a reduction of the width difference. 
The analyses of measurements of the width difference using B s decays into 
CP eigenstates have to be modified in the presence of new physics. We discuss 
how the width difference and the new CP violating phase in B s mixing can 
be measured. 
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I. INTRODUCTION 



The physics of the B s system is entering an exciting era, when experiments start exploring 
the range of mixing that is relevant to the Standard Model (SM) or to models of New Physics 
(NP). Recently, Dunietz (jT) has studied the possibility of using untagged B s samples for 
various measurements. The important ingredient in such an analysis is the large width 
difference that is expected in the B s system. In some sense, the B s system is similar to the 
neutral kaon system, where experiments study mass eigenstates {Kl and Kg) , rather than 
the Bd system, where experiments study the flavor eigenstates (B° and B°). 

In [|TJ, several ideas were put forward on how to measure the width difference, Ar, and 
how to extract the CP violating phase 7. However, the entire analysis was carried out 
in the framework of the SM. In this paper we investigate the implications of NP which 
can significantly affect the B s mass difference, AM. Phenomenological constraints on the 
relevant NP are rather weak. Moreover, various theories of flavor physics suggest that NP 
effects are more likely to appear in processes involving the heavy generation; the B s system 
is unique in that it is the only neutral meson that is expected to exhibit mixing and does 
not involve first generation quarks. 

We focus our attention on NP that contributes only to AM. This is the case in most 
extensions of the SM: new, heavy particles cannot compete with the iy-mediated tree level 
B s decay amplitudes, and thus to a very good approximation all decay amplitudes are given 
by the SM. However, the new effects could even dominate over the highly suppressed box 
diagram mixing amplitude. Naively, one would think that since in such models the B s decays 
are described by the SM, also the width difference is given by the SM. As we will discuss in 
this paper, this is incorrect. 

The organization of the paper is as follows: in Section II we collect the relevant formalism. 
In Section III we discuss how new sources of CP violation that in general arise in NP models 
reduce Ar. In Section IV we describe how the width difference can be determined in a 
model independent way, and we analyze several methods to extract the new CP violating 
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phase using untagged B s samples. We also argue that the unitarity angle 7 can be extracted 
even in the presence of NP. Finally, Section V contains our conclusions. 



II. FORMALISM 

We start by collecting the relevant formulae and definitions. Most of them are well known 
||. An arbitrary neutral S s -meson state, 



a\B s ) + b\ B.) 



(2.1) 



is governed by the time-dependent Schrodinger equation 
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Here M and T are 2x2 Hermitian matrices. CPT invariance guarantees H\\ 
eigenstates of the mass matrix are 



(2.2) 



H 22 . The 



Bl) = V I B 8 ) + q I B 8 ) , I B H ) = p \ B s ) - q \ B s ) , 



(2.3) 



with eigenvalues 



(al,h — Ml,h — - ^l,h ■ 



(2.4) 



Here Ml,h and denote the masses and decay widths of Bl^h and, by definition, Mr > 
Ml. In the SM > Th and then Bl is the short lived mass eigenstate. However, in the 
presence of NP this is not necessarily true. Therefore, we also label the mass eigenstates by 
B iong and B short , where, by definition, T short > Ti ong . We define 



AM = Mjj — Ml , Ah I' 



short L long j 



(2.5) 



The coefficient ratio q/p satisfies 



' q \ 2 _ (M* 2 -|ry 

v " (M 12 -|r 12 ) 



(2.6) 



The width difference is 

4|fle(M 12 r; 2 )| 

AM v ; 

The experimental lower bound AM/ T > 8.8 ||, implies AM ^> Ar and consequently |2j 
|Mi 2 | ^> | Ti2 1 ■ Thus, to a very good approximation 

AM = 2|M 12 |, |g/p| = l. (2.8) 

For a final state / we define the interference terms 

We further assume |(/ | S s )| = |(/ | S s )| (namely, no CP violation in decay) and then 
|A| = |A|. We always choose / such that |A| < 1. The rapid time-dependent oscillations, 
which depend on AMt, cancel in untagged data samples |]. We define 

r \f (t)} = r(B phys (t) -> /) + r(B phys (t) - /) , (2.10) 

where B p h ys (B p h ys ) denotes a time-evolved initially unmixed B s (B s ). We have 

r [/(*)] = r (^^/) |(i + |A| 2 ) (e^ + e- r ^)+2i?eA (e"™ - e" r ^) ) , (2.11) 
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J(t)] = r( ^ 2 "" /} |(l + |A| 2 ) (e-^ t + e- r « t )+2Re\ ( 



where we used |A| = |A| and T(B S — > /) = r(i? s — > /). Two limits are of particular interest. 
For decay modes which are specific of B s or B s , A = A = (semileptonic decays are flavor 
specific; b — ► cud decays are also flavor specific to a very good approximation), 

T [/ (f)] = r [J = r( ^ "" /} je"™ + e~ r -4 . (2.12) 



For a final state that is a CP eigenstate, A* = A with |A| = 1, and 

r (b s -> /) 



r [/ (t)] = 2 {e" rLt + e~ r «* + i?eA (e" r ^ - e"^') } . (2.13) 
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III. Ar IN THE SM AND BEYOND 



Ar is produced by decay channels that are common to B s and B s . In the SM, Ar 
is relatively large since the CKM unsuppressed parton decay b — > ccs produces such final 
states. Two approaches have been used to actually estimate Ar. The first is a quark level 
calculation, where Ar corresponds to the imaginary part of the box diagram. Assuming 
factorization ||, and including QCD corrections 0, the result is 

In deriving (|3.1| ), Ar = 2r 12 has been used, which is justified when Cr violation is neglected. 



This is a very good approximation in the SM. The second approach sums over exclusive 
channels ||. Assuming factorization and Heavy Quark Symmetry relations between form 
factors, 

Ar , . 

— «0.15. (3.2) 

The explicit dependence on the decay constants was not given in ||. In general the width 
difference scales like the decay constant squared. Therefore, the use of the recently measured 
fo s = 273 ± 39MeV J7|, instead of fo s = 230 MeV as used in || increases the estimate of 

Ar. 

When Nr contributes comparably to or dominates over the SM contribution to the B s 
mixing, Cr may be significantly violated. We now show that new Cr violating contributions 
to the mixing always reduce Ar relative to the SM prediction. Eqs. fl2.7|) and fl2.8|) lead to 



Ar = 2|r 12 | | cos 2^| , 2£ = arg(-M 12 r* 2 ) . (3.3) 

Under the reasonable assumption that Nr does not significantly affect the leading decay 
processes, Ti 2 arises from T(b — > ccs). Consequently, 2£ is the phase difference between the 
total mixing amplitude and the b — > ccs decay amplitude. In the SM, 

^'H-t§H (3 - 4) 



and then cos 2^ = 1 to a very high accuracy. With NP, new phases could be present, leading 
to cos 2^ < 1. This proves our statement. 

The reduction of Ar can be understood intuitively as follows. In the absence of CP 
violation, the two mass eigenstates are also CP eigenstates. The large Ar is an indication 
that most of the b — > ccs decays are into CP even final states. With CP violation, in the 
basis where the b — > ccs amplitude is real, the mass eigenstates are no longer approximate 
CP eigenstates. Then, both mass eigenstates decay into CP even final states. Consequently, 
Ar is reduced. 

The NP effects on the mixing amplitude M\% can be parameterized as 

M 12 = Mf 2 M (l + ac w ) , (3.5) 

so that £ is modified to (we approximate (3' = 0) 

/ 1 ( asin# \ 

2£ = arg (l + ae ie ) = tan" 1 — . 3.6 

v 7 \ 1 + a cos 9 ) 

The parameters a and 9 give the relative magnitude and relative phase of the NP contri- 
bution, i.e. a = \M^ 2 P /M^ M | and 9 = arg (^M^ p /M^^j. Phenomenological constraints do 
not exclude a > 1 [§], and 9 can have any value. Consequently, in the presence of NP, £ 
could assume any value. 

We do not study any particular NP model in detail, but mention here a few cases which 
predict (or can accommodate) large a and 9. In models with vectorlike down type quarks, 
a tree level bsZ vertex is induced. The experimental bounds on BR(B X s fi + fi~) || and 
BR(B — > X s vv) |H| allow for a < 0.25 and arbitrary 9 ||11|| . In fourth generation models 
the t' box diagram can be large and a > 1 with arbitrary 9 is allowed [j3|,|l2"|. In SUSY 
models without R parity, tree level sneutrino exchange (induced by the AL = 1 couplings) 
contributes to B s mixing. Then, a > 1 with arbitrary 9 is allowed |p~3f . 

Finally, we mention the possibility of NP that affects the decay rates. An example would 
be models that give large BR{b — > sg) = O(10)% [TJJ]. Since b — ► sg can produce CP 
eigenstates, Ti2 could change as well. In particular, Ar may become somewhat larger than 
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the SM prediction. For the rest of the paper, we restrict ourselves to the more reasonable 
case, namely when only M 12 is significantly modified by the NP. 

To conclude this section: If NP contributes significantly to the B s mixing, the likely 
results are the enhancement of AM and a suppression of AI\ If experiments are able to 
push the lower bound on AM up, or the upper bound on Ar down, they might, in principle, 
be able to find NP. 

IV. MEASURING Ar, f AND 7 

In this section, we examine how to measure Ar and the angles £ and 7 using untagged 
B s data samples. In Ref. [|J, several methods of how to measure Ar and 7 have been 
analyzed in the framework of the SM. Here, we study how they have to be modified in the 
presence of NP and how they can probe it. We start by looking to the methods that were 
proposed to extract Ar. We show that by combining few measurements, both Ar and £ 
can be extracted. 

The first method is to fit the time dependent rate of flavor specific decays (see Eq. ( |2.12| )) 
to two exponentials. This fit determines Ti mg and T s h ort) and therefore, V and Ar. This 
method is still valid in the presence of NP. 

The second class of methods to measure Ar uses decays into CP eigenstates. (Unless 
otherwise specified, when we talk about CP eigenstates we refer to those that are produced 
by the b — > ccs decay.) For concreteness, we focus on one method. We assume that T is 
known from either a fit to a flavor specific data samples, or from the average b hadrons 
lifetime IJ. Then, in the SM, the B s lifetime measured using B s decays into CP even final 
states determines T s h or t, and therefore, Ar. (Examples of such final states are D+Dj and 
the CP even component of J/ip<f) that can be obtained using transversity analysis [|15j.) As 
mentioned in these methods are modified in the presence of NP. We now study this point 
in detail. 

In a general NP model, when decays into CP eigenstates are dominated by the W 
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exchange tree diagram, we have 

A = ±e 2i «, (4.1) 

where, by convention, +(— ) stand for a CP even (odd) final state. We emphasize that the 
phase £ in (|4.1|) and ( |3.3p is the same phase. Then the decay rate into CP eigenstates fl2.13| ) 
satisfies 



T(B -> CP even, t) oc cos 2 £ e~ l + sin^ e" 1 H \ (4.2) 

For a decay into a CP odd state, and Th are interchanged. 

A comment about discrete ambiguities is in order. In the SM the light state has a shorter 
lifetime and it is a CP even state |[|. However, in the presence of NP things may be different. 
When the NP contribution is large and negative (cos 2£ < 0) the heavy state has a shorter 
lifetime, thus reversing the SM prediction. For sin 2£ ^ the mass eigenstates are no longer 
CP eigenstates, both eigenstates decay into CP even and CP odd states (see Eq. (|4.2|) ). 
We emphasize, however, that always the longer lived state is "closer" to the CP odd state, 
namely, its decay width into CP odd states is larger than the decay width of the other mass 
eigenstate. Therefore, because AM/M is tiny, sign(r i — T H ) = sign(cos2£) is practically 
undetectable. Moreover, because 2£ always appears as an argument of a cosine function, 
also sign(sin2£) cannot be determined using untagged B s samples. (Note, however, that if 
the rapid AMt time dependent CP asymmetries in B s decays will be measured, this sign 
will be determined, since asymmetries in modes governed by b — > ccs can be used to extract 
sin2£ 0.) In conclusion, in our analysis, there is a four fold ambiguity in the value of 2£. 



A. Combining flavor specific and CP eigenstates data 

In principle, a three parameter fit of a decay into a CP even eigenstate can be used to 
measure T, Ar and £ using Eq. ( (4.2Q . Even if this cannot be done in practice, by comparing 
the measurements of Ar from flavor specific decays and CP eigenstate decays, £ can be 
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measured. Experimentally, most of the data are expected to be taken for small T t. Then, 
using Aft < 1, Eq. (O) becomes 



r^CPeven^oce-^ 4 , r + = ( T + ^£i^Mj . (4.3) 

Using T and Ar as measured from the flavor specific data, a one parameter fit to the decay 
rate gives £. Actually, such a fit determines 

Ar CP = 2(r+ - r) = Ar 1 cos 2f | . (4.4) 

By comparing it to the real width difference as measured using flavor specific data, ATps, 
we get 

AF 

|cos2£| = — ££. (4.5) 

Al FS 

This method would be particularly useful if £ is neither very small nor very large. For 



£ ~ 7r/4 the width difference becomes too small to be measured (see Eq. (|3.3|) ). For £ ~ 
the precision of the measurement should be very high. While we concentrate on one example, 
we emphasize that Eq. (|4.5| ) is relevant to all the methods for measuring Ar using decays 
into CP eigenstates as suggested in 



B. Theory 

Since £ = in the SM, we get from Eq. (|3.3|) 

Ar = Ar SA/ |cos2f|. (4.6) 

We learn that, if we knew the SM prediction for AT$m, the measurement of Ar would allow 
for the determination of £. The problem is, of course, that we do not know how to calculate 
ArsM to a high accuracy. The box diagram calculation suffers from the uncertainty in 
fs a and, in any case, can be used only as an order of magnitude estimate since hadronic 
physics alters it. Similarly, there are large uncertainties in using the sum over exclusive 
states. However, we are able to get bounds on £ even with these large uncertainties. To 
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do that, we need only rough bounds on the SM range. The upper bound is rather safe, as 
Ar < F(b — ► ccs). Based on the inclusive ]TB| and the exclusive calculations, we can 



conservatively take an upper bound of BR(b — ► ccs) < 40%. The lower bound is less reliable. 
Even if one can obtain a lower bound on BR(b — > ccs) it would not help since we only know 
that Ar/r < BR(b — > ccs). However, from Eqs. ( |3.1p and ( |3.2| ) we can conservatively take 
Ar/r > 10%. We believe then that it is safe to assume 

Ar 

10% < < 40% . (4.7) 

Therefore, a measurement of Ar will provided an upper bound on £: 

and similarly, any upper bound on Ar below the SM prediction, would give a lower bound 
on £: 

cos2?< (^Ll5%' (4 ' 9) 

Note that Ar/r < 0.1 would be evidence for NP; this is an unusual situation. Usually, NP 
enhances observables. In this case the SM is at the maximum and NP can only reduce Ar. 

C. Time tag 

It will be interesting if very late decaying B s samples can be collected [|TJ. Then, all 
the short lived B s have already decayed and the beam consists purely of the long lived 
state, Bi on g. Of course, for any useful measurement, large statistics is needed. The expected 
number of reconstructed B s — > J/ip<p events is about 1.2 x 10 4 at CDF in run II ||18|| , and 
about 3.8 x 10 5 per year at the LHC |T!|. Assuming that this can be achieved, £ could be 
measured using the time tag method. 

From Eqs. fl2TT2D and ( gig ) we get 



R (t) = oc (i _ Re X tanh ( —) \ , ( 1.10) 



r[g(t)\ I V 2 
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where / is a CP even eigenstate and g is a flavor specific state. To enhance statistics, a sum 
over several CP even final states can be performed. Of course, if the time dependence can 
be traced, Ar and £ can be extracted (recall, ReX = cos2£). In practice it might be easier 
just to make a cut to select events at large t, when tanh (Art/2) — > 1. This can be thought 
of as an almost pure Bi ong beam. Then £ can be obtained from the ratio 

Finally, let us estimate the loss in statistics if this method is used. We define 

P{t) = — - N{Blon9 ;j) n r = (l + exp(-Art)r\ (4.12) 



N(B short ,t) + N(B long ,t) 

where N(B a ,t) oc exp(— T a t) is the number of B a in the beam (a = long, short). In order 
to get a data sample with purity P for a given Ar we need to make the cut at t cut such that 



r t r ,,t > 



The loss in statistics is exp(— r t cut ). For example, for P = 92% and AT/T = 0.30 one has to 
wait about 8 lifetimes, a loss of about 3 x 10 3 in statistics. We learn that it may be possible 
to use the time tag method. 

D. Measuring 7 

In |l[ several methods had been proposed for extracting 

7 = arg f_M^y (4.i4) 

Here we briefly describe the methods and show how they are modified in the presence of 
NP. More details on the methods, including the validity of the assumptions that are used, 
are given in Ref. 

There are basically two classes of methods. The first one is based on decays into CP 
eigenstates mediated by the b uud transition, e.g. p°Ks- When T and Ar are known, a 
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one parameter fit to the time dependence of untagged p°K$ events determines ReX (see Eq. 
( P-13| )). Assuming that the penguin contamination is small, 

q (p°K s | B a ) 

p(p°K s \B s ) ■ l4 - i5j 

Then, ReX = — cos(2£ + 27) and £ + 7 can be extracted. 

The second class uses pairs of final states that are CP conjugate of each other and are 
mediated by the b —* cus and b — > ucs transitions, e.g. D~K + and DfK~. When the total 
decay width, T(B S — > f) is known a three parameter fit to the time dependences of the 
decays into / and / determines |A|, ReX and ReX (see Eqs. fl2.11|) ). In general, 

A -p(D 7 K+\B s ) |A|e ' X ~ q{DfK-\B s )-^ e ' ^ ibj 

where 5 is the strong phase. Then, ReX = \X\ cos(<5 + 2^ + 7), ReX = \X \ cos(5 — 2^ — 7) and 
2^ + 7 can be extracted. 

We learn that when T, Ar and ^ are known (from the combination of the flavor specific 
and the b ccs CP eigenstates decays data) all these methods can be used to measure 7 
up to discrete ambiguities. 



V. CONCLUSIONS 

Various extensions of the Standard Model predict new significant contributions to the 
B s mass difference. In general, such models introduce new sources of CP violation beyond 
the single phase of the CKM matrix. Still, in most cases, the B s decay rates are described 
by the SM. Nevertheless, the width difference can be significantly reduced. Actually, such a 
reduction is an indication of CP violation: the large SM prediction for Ar is based on the 
fact that the decay width into CP even final states is larger than into CP odd final states. 
When new CP violating phases appear in the mixing amplitude then the mass eigenstates 
can sizably differ from the CP eigenstates, and both mass eigenstates are allowed to decay 
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into CP even final statesf]. Consequently, Ar reduced. 

Since in the SM the CP violating phase £ is very small 0, (£ < |T4sKi&AcsK*| < 
2.5 x 10~ 2 ), a measurement of £ 7^ would be a clear evidence for CP violation beyond the 
SM. Actually, £ 7^ is the equivalent of a + j3 + j 7^ tt, which seems to be much harder to test 
experimentally |2U|. Furthermore, £ can be extracted using the B s leading b — > ccs decay 
modes (CKM unsuppressed) . All these reasons explain why the measurement of £ is very 
important. If the rapid AMt oscillation can be traced, the time dependent CP asymmetries 
in B s decay modes mediated by b — > ccs will measure sin 2£. However, untagged B s samples 
seem to be the best place for this measurement. We described several ways of determine £ 
using untagged B s samples: 

• Compare the SM calculation (which has large uncertainties) to the measured Ar, 

|cos2£| = — — . (5.1) 

^1 SM 

• Compare Ar as measured from decays into CP eigenstates, to Ar from flavor specific 
decays, 

|cos2£| = ^. (5.2) 

• If a measurement at late times is possible, decay rates of a specific mass eigenstate 
can be measured, and then 

os rW^CPeve^)^ 
1 [B s — > CP even) 



*When \q/p\ = 1 one can always choose a convention for CP transformation such that the mass 
eigenstates are CP eigenstates. But then, CP is violated by the decay amplitudes and both CP 
eigenstates decay into CP even states. Of course, the final result is the same, no matter what 
convention we use. 
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